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Abstract. We investigate the angular correlation between pions and photons produced in deuteron-gold 
collisions at RHIC and proton-lead collisions at LHC using the Color Glass Condensate formalism and 
make predictions for the dependence of the production cross section on the angle between the pion and 
photon at different rapidities and transverse momenta. Measuring this dependence would shed further light 
on the role of high gluon density and saturation dynamics at RHIC and LHC. 
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1 Introduction 

The Color Glass Condensate (CGC) formalism has been 
quite successful in describing single inclusive production as 
well as multiplicity of hadrons at RHIC [lj. While CGC 
based predictions for the suppression of forward rapid- 
ity spectrum of produced hadrons in deuteron-gold col- 
lisions were observed in the experimental data [2], there 
have been more recent and more traditional approaches 
which can now also fit the data [3]. Therefore, it is desir- 
able to investigate more processes which could help ver- 
ify/falsify the claim that saturation dynamics has been 
conclusively discovered at RHIC . Electromagnetic pro- 
cesses [4] are specially promising tools in investigating sat- 
uration physics since they do not suffer from the theoret- 
ical uncertainties associated with hadronization. Another 
venue which is particularly sensitive to saturation effects, 
and can serve a powerful discriminator between different 
approaches is two particle correlations [5]. 

Two particle correlations are expected to be a sensitive 
probe of high energy QCD and saturation dynamics where 
one expects a decorrclation between the two hadrons when 
the rapidity separation between the two hadrons increases. 
This is due to the large phase space (in x or rapidity) 
between the produced hadrons. Since the small x evolu- 
tion involves a random walk in transverse momentum, this 
washes away the angular correlation which is the opposite 
of the case in pQCD in the DGLAP limit where the two 
produced hadrons are back to back. 

While two hadron production cross section has been in- 
vestigated using the CGC formalism and analytic results 
are known [6], the resulting expressions are complicated, 
specially when the two hadrons are separated by large 
rapidities. In order to get results which lend themselves 
easily to numerical methods, one needs to make further 
assumptions and approximations which reduce the accu- 
racy of the results and diminish the predictability of the 
formalism. For instance, the two hadron production cross 
section involves higher point functions of Wilson lines un- 
like the single inclusive production which involves only the 



two point function. One would then need to either solve 
the complicated JIMWLK PQ equations for these higher 
point functions or make a mean field (or Gaussian) ap- 
proximation where all higher point functions are approx- 
imated as powers of the two point function. This would 
clearly neglect some correlations which may be important. 

Here we investigate photon-pion production in deuteron- 
gold (or proton-lead) collisions which is a much simpler 
problem than two hadron production but still has the es- 
sential physics of saturation. The basic building block in 
this process is the quark anti-quark dipole scattering prob- 
ability, denoted Np, which also shows up in single inclu- 
sive photon or dilepton production, as well as single in- 
clusive hadron production in the forward rapidity region. 
The dipole scattering probability satisfies the JIMWLK 
(or BK in large N c and mean field approximation) and its 
properties are understood quite well even. Even though an 
analytic solution to the dipole evolution equations is not 
known, it can be solved numerically 7 and approximate 
analytic solutions are known in various limit. 



2 The production cross section 

The cross section for inclusive production of a photon and 
a parton in the scattering of a quark on a given target T 
such as a proton or nucleus is given by [4] 
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where It, k t are the transverse momenta of the final state 
quark and photon respectively, z is the fraction of the 
energy of the initial state quark carried by the photon and 
x g is the momentum fraction of the gluons in the target 
T (proton or nucleus) given by 
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where y q and y 7 are the rapidities of the final state quark 
and photon. To relate (p} to the production cross section 
for a hadron and a photon in a deuteron (proton)-target 
scattering, we will convolute |T]) with quark distribution 
functions of a proton (or deuteron, using isospin symme- 
try) given by CTEQ6 and pion fragmentation functions 
given by KKP and rewrite fl} as (below we will consider 
impact parameter integrated cross sections only) 
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with qt — zlt and where a 



is a angle 



2 v / 2(2tt) 3 V« 

independent constant which will drop out when we con- 
sider the angular correlations below and qt, q~ are the 
transverse momentum and energy of the produced pion. 
The angle 9 is defined to be the angle between the pro- 
duced photon and hadron. The lower limit of the z in- 
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tegration is given by z r , 
x p = -±=(q t /z e~ yh +k t e" 

energies and energies of the produced pion are related via 
q~ = e Vh and similarly for the photon. 

Before a quantitative numerical analysis it is worth 
looking at some limits of (JT]) . It can be shown [J] that it the 
limit where the photon has a very large transverse momen- 
tum (such that the collinear photon-quark contributions 
are suppressed), one recovers the standard pQCD result 
for q g — ► q 7 convoluted with the gluon distribution func- 
tion of the target. On the other hand and in the collinear 
limit, one recovers the fragmentation photons [8]. There- 
fore, the above expressions include contributions of both 
direct and fragmentation photons. On the other hand, it 
does not include the contribution of gluons from the pro- 
jectile scattering from the target and then fragmenting 
into photons which necessarily involves quark production 
in the intermediate stage and is therefore higher order in 
a s even though it may be numerically important in mid 
rapidity region. For the sake of theoretical consistency, 
this contribution is not included here. 

Therefore, for the case when the produced pion and 
photon have the same transverse momentum and rapidity, 
one expects a collinear divergence as the angle 9^0 
in the correlation function which disappears as one goes 
to larger angles where in the limit where 9 — > 180 one 
recovers the standard back to back configuration. If the 
pion and photon are separated by a large rapidity while 
at the same transverse momentum, the cross section does 
not receive any contribution from collinear configurations. 



2.1 Numerical results 

In order to investigate the angular correlation between 
the produced pion and photon, we define the correlation 



function C{9) as 
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We expect that C(9) — > 1 when 9 — > when the pion and 
photon have the same rapidity and transverse momentum 
since the cross section in this case is dominated by the 
collinear divergence between the photon and pion. 

In Fig. ([!} we show the correlation function C(9) when 
both photon and pion are in mid rapidity and for several 
transverse momenta, for deuteron-gold collisions at RHIC. 
It goes to 1 when the angle 9 goes to while at 9 — > 180, 
its magnitude increases as one goes to higher transverse 
momentum. We note that the magnitude of the correlation 
function changes by almost an order of magnitude as one 
increases the transverse momenta of the pion and photon 
from 2 GeV to 10 GeV when 9 — * 180 where the qg — + (77 
process dominates. This increase of the correlation func- 
tion is due to the fact that as one goes to higher transverse 
momenta, high gluon density effects become weaker. 




Fig. 1. Pion- photon correlation function in deuteron-gold col- 
lisions at RHIC. 



In Fig. ([2|) we show the correlation function C(9) in 
the case when the photon and pion are separated by a 
large rapidity (yh — 4, y 1 = 0) but are at the same trans- 
verse momenta. As discussed earlier, the contribution of 
fragmentation photons becomes negligible and the correla- 
tion function is dominated by the quark-gluon scattering. 
Unfortunately, the transverse momentum reach in the for- 
ward rapidity region of RHIC is limited to low transverse 
momenta due to the kinematics. We also note that the 
small dip in the region when 9 —* 180 is numerical artifact 
and is caused by the rapid oscillation of the dipole profile 
which makes it extremely difficult to Fourier transform it 
accurately Q. In Fig. ([3|) we show the correlation function 
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when the two particles are separated by a large rapidity 
but with yh — and y 1 = 4. We also show the effect of an 
isolation cut on the photon where the cut parameter e is 
independent of the cone size at Leading Order and is taken 
to be e = 0.1 at RHIC. This cut affects the minimum of 
the z integration in the cross section to 



(5) 



It is clear that the isolation cut would severely diminish 
the cross section unless the rapidity separation between 
the pion and photon is large and negative, or when the 
transverse momenta of the two particles differ by orders 
of magnitude. Therefore, the isolation cut is not applied 
to the results shown in the other figures. 



Finally in Fig. (0} we show the correlation function for 
proton- lead collisions at LHC energies (y/s — 8800 GeV). 
A similar pattern in the angle dependence is observed ex- 
cept that the transverse momentum dependence of the 
correlation function is much weaker in the back to back 
case. This is presumably due to the larger kinematic re- 
gion in LHC where leading twist shadowing is more impor- 
tant as compared to RHIC. In summary, hadron-photon 
angular correlation function provides an additional tool 
with which to investigate saturation dynamics at RHIC 
and LHC. It depends on the same basic ingredient as pho- 
ton and dilepton production and single hadron production 
in the forward rapidity region. Therefore an experimen- 
tal verification/falsification of this correlation would shed 
much light on the domain of applicability of CGC and 
saturation physics. 
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Fig. 2. Same as in Fig. |T} but with different rapidities. 
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Fig. 4. Pion-photon correlation function in mid rapidity 
proton-lead collisions at LHC. 
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Fig. 3. Same as in Fig. ([T]) but with different rapidities and 
with or without an isolation cut. 
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